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A static scission configuration in cold ternary fission has been considered in the 
framework of two mean field approaches. The virial theorems has been suggested 
to investigate correlations in the phase space, starting from a kinetic equation. The 
inverse mean field method is applied to solve single-particle Schrodinger equation, 
instead of constrained sclfconsistent Hartree— Fock equations. It is shown, that it 
is possible to simulate one-dimensional three-center system via inverse scattering 
method in the approximation of reflectless single-particle potentials. 

1 Motivation 

Ternary fission involving the emission of a-particle was first observed more 
than fifty years ago H. Emission of a-particles in the spontaneous fission of 
™ Cf has also a long history of investigation experimentally a and theoretically 
uu. A renewed interest in these processes arose in connection with 

• modern experimental techniques (7—7—7 and x— 7— 7 triple coincidence, 
(Gammasphere with 110 Compton suppressed Ge detectors), which allow 
the fine resolution- of the mass, charge and angular momentum content 
of the fragments. p 

• Recently direct experimental evidence was presented for the cold (neu- 
tronless) ternary spontaneous fission of 252 Cf in which the third particle 
is an a-particle la, or 10 Bel3. This confirms that a large variety of nuclear 
large-amnlit|Ude-|Collective motions such as bipoeLal fission @, cold binary 
fission BEli^'Oli 2 ], heavy cluster radioactivity t2rli-3, and inverse processes, 
such as subbarier fusion 113, could belong to the general phenomenon of 
cold nuclear fragmentation. 

• Cluster like models were used successfully to reproduce general fea- 
tures of the cold ternary fragmentation. However the scission configura- 
tion has been built in fact by hands. 

Therefore it is actual to develop microscopical or semi-microscopical ap- 
proach to this scission-point concept of nuclear fragmentation. There are 
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well developed methods to calculate, in the framework of many-body self- 
consistent approach, static properties of a well isolated nucleus iu its ground 
state. There also exists a well developed two-center shell model However, 
a three^center shell model has not been developed yet, except for very early 
steps Three-center shapes are practically not investigated, in compari- 
son with the two-center ones. There-exists the generalizations of mean-field 
models to the case of two-centers E3, but a ternary configuration is out of 
consideration, because of uncertainties to select a peculiar set of constraints. 

There exist a number of calculations for nucleus-nucleus collisions in the 
frame of time-dependent mean-field methods, but an evolution of the cold 
fragmentation has not been investigated yet. Therefore, although the prin- 
cipal way to describe nuclear fragmentation in the framework of many-body 
self-consistent approach exists, it is interesting to develop other mean-field 
approaches to analyse these phenomena from different points of view. 

In this Contribution we suggest two methods to analyse a static scission 
configuration in cold ternary fission in the framework of mean field approach. 
In Section 2, starting from kinetic equation the virial theorems has been sug- 
gested to investigate correlations in the phase space. It gives possibility to 
formulate in future generalised set of constraints in momentum and cartesian 
spaces. In Section 3 the inverse mean field method is applied to solve single- 
particle Schrodinger equation, instead of constrained selfconsistent Hartree- 
Fock equations. It makes irt possible to simulate one-dimensional three-center 
system via inverse scattering method. 

2 Virial theorems 

Within the mean-field approximation, we analyze the evolution of one-body 
Wigner phase-space distribution function of the full many-body wave function, 
following the well developed scheme using as the starting-point the Vlasov 
equation for the Wigner phase-space distribution function c3 

dj_ + I d£ _ dV dj_ = 

dt m df df dp re ' 

with a "relaxation term" I re i is added to the kinetic equation to describe 
dissipation effects. The quantity V(r, t) is the self-consistent single-particle 
potential which is assumed here to be local, m is the mass of nucleon. 
Out of the Wigner distribution function, virials at different orders can help 
to extract useful physical information from the total phase space dynamics. 
Integrating the initial kinetic equation (|l|) over the momentum space with 
different polynomial weighting functions of the p-variable one comes, as well 
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know E3,EjjHil, to an infinite chain of equations for local collective observables 
including the density, collective velocity, pressure and an infinite set of tenso- 
rial functions of the time and space coordinates, which arc defined as moments 
of the distribution function in the momentum space: 

• the particle n(f,t) = J dp f (r, p, t) , and the mass p(r,t) = mn(r, t) den- 
sities, 

• the collective current and velocity of nuclear matter 

p(r,t)u(r,t) = / dppf{f,p,t), 



• the pressure tensor and the energy and momentum transfer tensors of 
different orders 

¥ij(r, t) = — / dpq^f^p, t), q t =p % - mUi, 
m J 

P ij..k(^ t ) = —^T / d P mj-1kf(r,p,t), 



and the integrals related to relaxation terms 

dpl rei =0, / dppl re i = 0, 



lij = — I dpqiq^Irei, 



III 



Truncating this chain at order two in q one arrives at the "fluid dynamical" 
level of description of nuclear processes. 

k 

Du % x - d¥ ik v 2 

k k 

+ p^2 £isj{2Q s Uj + -jf x j) = °) ( 3 ) 
Dt ^ V dxk oxk oxk 
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m 

dx k J \ dt 



— [ dp q l q ] I r eU 
m J 



where the usual notation -Sr = -|j + X) M fc d§~~ ^ s introduced for the operator 

giving the material derivative, or the rate of change at a point moving locally 

with the fluid. The hydrodynamical set of Eqs. (|||5|) describes an evolution of 

a rotating nuclear system. We consider two frames of reference with a common 

origin: an inertial frame, (Xi, X2, X3), and a moving frame, {x\, X2,Xs). Let 
3 

Xi = "^ijXj be the linear transformation that relates the coordinates, X 

and x, of a point in two frames. The orientation of the moving frame, with 
respect to the inertial frame, will be assumed to be time dependent. Since 
Tij (t) must represent an orthogonal transformation, the vector 



1 j,k,m V m J )'» 

represents a general time-dependent rotation of the af-frame with respect to 
the inertial frame. 

Let us define integral collective " observables" (the integrals over the whole 
phase space of one nucleon containing the distribution function appropriatly 
weighted), namely an inertia tensor §ij(t), the dynamical part of the angular 
momentum Li (t) , the integral pressure tensor Iljj (t) defined as 



dxXiXjp, Uij — J (ixPjj, 

Lk = / dx pxiUj, 

The dynamics in terms of the latter ilpbservables" is expressed by a set of 
virial equations in the rotating frame cB 



d 2 



13 



2^2^ s J df pu k {e isk Xj + EjskXi) 



extr2000: submitted to World Scientific on February 8, 2008 4 



2W.y - ZKij - 2IIy 



—jr (eisfeJfcj + £jsfcJfci) = 0, 



dt 



pUkX s 



i,j,m s 

— Ily + Fy + 2 ^ fi s (e j s fcllfc.,- + ejsfcllfcj) = Ry, 

where the tensors of collective kinetic and potential energies, and the relax- 
ation tensor are 

f _ f ^ dV 

Ky = dx UiUjp, Wij = J dx Xj-g—n, 

The above equations constitute a formal framework within which the coupling 
of the deformations in the r-space and in the p-space can be explicitly worked 
out. It may help to fofmulate in a future a generalised set of constrains in the 
phase space. 



3 Inverse mean field method 

3.1 The framework 

Methods of nonlinear dynamics, gave yet the possibility to derive for nuclear 
physics unexpected collective modes, which can not be obtained by tradi- 
tional mr;thods_pf perturbation theory near some equilibrium state (see e.g. 
review E3 and c3 for the recent refs.). The most important reason is that the 
fragmentation and clusterization is a very general phenomenon. There are 
cluster objects in subnuclear and macro physics. Very different theoretical 
methods were developed in these fields. However, there are only few basic 
physical ideas, and most of the methods deal with nonlinear partial differ- 
ential equations. One of the PpQst important part of soliton theory is the 
inverse scattering method l3E3E!l and its applications to the integration of 
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nonlinear partial differential equations c3. The inverse methods to integrate 
nonlinear evolution equations are often more effective than a direct numerical 
integration. Let us demonstrate this statement for the following sinmle case. 
The type of systems under consideration are slabs of nuclear matter EJ, which 
are finite in the z coordinate and infinite and homogeneous in two transverse 
directions. The wave function for the slab geomethry is 

1 H 2 k 2 
^kj.n(x) = —=ip n (z) exp(ikj_r), e k±n = —± + e n , (5) 



where r = (x, y), \i± = (k x , k y ), and f2 is the transverse normalization area. 
h 2 d 2 

- ^^2^"( z ) + U(z)ip n (z) = e n ip n (z), (6) 

A direct method to solve the single-particle problem (||) is to assign a func- 
tional of interaction £ (usually an effective density dependent Skyrme force), 
to derive the ansatz for the one-body potential, as the first variation of a func- 
tional of interaction in density U{z) = U[p{z)\ = 5£/8p. Then to solve the 
Hartree-Fock problem under the set of constraints, which define the specifics 
of the nuclear system. In the simplest case of a ground state, one should 
conserve the total particle number of nucleons (A), which is related to the 
"thickness" of a slab, via A =*> A = (6A P 2 N /Try/ 3 , which gives the same 
radius for a three-dimensional system and its one-dimensional analogue. As 
a result, one obtains the energies of the single particle states e„, their wave 
functions ip n (z), the density profile p(x) =>■ p(z) 



No N 2m 

p{z) = ^2 a n4>l(z), -A. = ^ a W) a n = —^(e F - e„), (7) 

n—l n=l 



and the corresponding single-particle potential. a n are the occupation num- 
bers, No is the number of occupied bound orbitals. The Fermy-energy e_p 
controls the conservation of the total number of nucleons. The energy (per 
nucleon) of a system is given by 

77 fc2 foo j I No ^oc 

{Y.^j J^fdz + lY. a2 r)+^j jW)]dz. (8) 



A 2mA 

n—l 



Finally, the set of formulas (|5|~|8|) completely defines the direct self-consistent 
problem. Following the inverse scattering method, one reduces the main dif- 
ferential Schrodinger equation ([|) to the integral Gel'fand-Levitan-Marchenko 
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equation 

K(x,y)+B(x + y)+ I B{y + z)K(x, z)dz = 0. (9) 



for a function K(x, y) . The kernel B is determined by the reflection coefficients 
R(k)(ek = h 2 k 2 /2m), and by the N bound state eigenvalues 

n i r°° 

B(z) = C?(/c n ) + - R(k) e-np(ikz)dk, e„ = -h 2 n 2 /2m. 

n=l 

N is the total number of the bound orbitals. The coefficients C n are uniquely 
specified by the boundary conditions and the symmetry of the problem under 
consideration. The general solution, U(z) = — (fi 2 /m)(dK(z, z)/dz), should 
naturally contain both, contributions due to the continuum of the spectrum 
and to its discrete part. There seems to be no way to obtain the general 
solutionXI^z) in a closed form. Eqs. (q),(§) have to be solved only numerically. 
In Ref. E-3 we used the approximation of reflectless (R(k) = 0), symmetrical 
(U(—z) — U(z)) potentials. This gave the possibility to derive the following 
set of relations 

h 2 rP- 9fr 2 N 

U{Z) = " m dT 2 Mdet l|M|l) = £ 



n=l 



N 



^n{z) = y^(M l ) n ih{z), A„(z) = C n {n n ) exp(-K n z), 

n=l 

m t„\ x A„(z)A;(z) / ITT ^!L±_^ iV /2 /im 

M n l{z) — Onl H ; , Cn(Kn) = 2k„ . (10) 



Consequently, in the approximation of reflectless potentials (R(k) — 0), the 
wave functions, the single-particle potential and the density profiles are com- 
pletely defined by the bound state eigenvalues via formulas (P),(|l0|). 

3.2 Results and Discussion 

In Ref. a series of calculations for the different layers, imitating nuclear 
systems in their ground state was provided. For a direct part of the calcu- 
lations by the Hartree-Fock method, the interaction functional was chosen 
in the form of effective Skyrme forces. The calculated spectrum of bound 
states was fed into the scheme of the inverse scattering method, and the rela- 
tions were used to recover the wave functions of the states, the single-particle 
potentials, and the corresponding densities. In this note, we generalize this 
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Figure 1. The density profiles of the light (A as 20, A ss 1.0) three-levels (AT = 3, JV = 3) 
model system calculated in the frame of inverse mean field method, a) the ground state 
(solid line); a fragmented two-center configuration (dotted and dashed lines); b) a ternary 
fragmentation of the system into three fragments. 



method to the case of fragmented configuration, trying to imitate two- and 
three-center nuclear systems. We use here only the inverse mean-field scheme 
( pp| ) . The details of the approach and systematic calculations of fragmented 
nuclear systems will be provided in a forthcoming publication. In Fig.l we 
present the results of the calculations of three-level (N = 3, Nq = 3) light 
(A ~ 20, A ~ 1.0) model system simulating the ground state (Fig. |l|(a) 
solid line), and fragmentation of the system into two fragments (Fig. |l](a), the 
dotted and dashed lines). In the same figure (Fig. |l|(b)) we present the frag- 
mentation of the system into three fragments (solid and dotted lines). One 
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can see that it is possible to simulate a one-dimensional three-center system 
via inverse scattering method. The following conclusions can be drawn. 

• The density profiles, calculated in the framework of inverse method, are 
practically identical to the results of calculation by SHF method. These 
results are valid for the ground state and for the system in the external 
potential field. 

• The global properties of single-particle potentials (the depth and an ef- 
fective radius) have been reproduced quite well, but the inverse method 
yields the quite strongly pronounced oscillations of the potential distri- 
butions within the inner region, and slightly different asymptotic tails 
of potential. In the framework of inverse scattering method, all bound 
states are taken into account in the calculation of the potential (|To| ) , but 
for the density distribution only the occupied states are taken into ac- 
count (see Eqs. Therefore, the slope of the tails of the potential and 
of the density distributions will we different. 

• We used, the approximation of reflectionless potentials, which gave us 
the possibility to obtain a simple closed set of relations (flCl), to calcu- 
late wave functions, density distributions and single particle potentials. 
The omitted reflection terms (R(k) = 0)are not important for the evalu- 
ation of the density distributions, due to the fact that only the deepest 
occupied states are used to evaluate density distribution (see Eq. (^)). 
The introduction of these reflection terms will lead to a smoothing of the 
oscillations in the inner part of the potential and to a correction of its 
asymptotic behaviour. 

• The presented method gives a tool to simulate the various sets of the 
static excited states of the system. This method could be useful to pre- 
pare in a simple way an initial state for the dynamical calculations in the 
frame of mean-field methods. 

4 Conclusions 

Recent experimental progress in the investigation of cold nuclear fragmen- 
tation has made the development of theoretical many-body methods highly 
desirable. Modern variants of self-consistent Hartree-Fock and relativistic 
mean-field models give the principal way to describe nuclear fragmentation in 
the framework of many-body self-consistent approach. However, the general- 
ization of these approaches to three-center case is not provided yet because of 
existing difficulties to select a suitable set of constraints. 
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We suggest two methods to analyse a static scission configuration in cold 
ternary fission in the framework of mean field approach. The virial theorems 
has been suggested to investigate correlations in the phase space, starting from 
a kinetic equation. The inverse mean field method is applied to solve single- 
particle Schrodinger equation, instead of constrained selfconsistent Hartree- 
Fock equations. It is shown, that it is possible to simulate one-dimensional 
three-center system via inverse scattering method in the approximation of 
reflectless single-particle potentials. 

These models may be useful as a guide to understand the general prop- 
erties of fragmented systems and to formulate the suitable set of constraints 
for the realistic three-dimensional mean field calculations of the three-center 
nuclear system. 
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